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Abstract 

In this paper, interpolation by scaled multi-integer translates of Gaussian kernels is studied. 
The main result establishes L p Sobolev error estimates and shows that the error is controlled 
by the L p multiplier norm of a Fourier multiplier closely related to the cardinal interpolant, and 
comparable to the Hilbert transform. Consequently, its multiplier norm is bounded independent 
of the grid spacing when 1 < p < oo, and involves a logarithmic term when p — 1 or oo. 



1 Introduction 

We consider interpolation by means of linear combinations of translates of a fixed Gaussian where 
the data consists of samples f(hk), k in Z n , of a continuous functions / and derive L p error estimates 
in terms of h and appropriate smoothness properties of /. Namely, our interpolants are of the form 

s h (x) = a k e-\ x ~ hk \ 2 

and, roughly speaking, our estimates are of the form ||/ — s^||r < C7i fe ||/|| W fc, valid for sufficiently 
large k. 

In a series of papers \18\ [T^l [i"5j 116] . Riemenschneider and Sivakumar have developed a compre- 
hensive theory of cardinal interpolation by Gaussians, treating issues of existence/uniqueness of 
interpolants, decay of fundamental functions, bounds on Lebesgue constants, and L p stability for 
data in £ p . This is an outgrowth of a general theory of cardinal interpolation that started for 
univariate splines with Schoenberg [T7], was extended to several variables for box splines by de 
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Boor, Hollig, and Riemenschneider [5], and for radial basis functions by Madych and Nelson [9J, 
and Buhmann[2]. 

One important topic not addressed by Riemenschneider and Sivakumar is the rate of convergence 
of the interpolant to a given smooth target function. In general, global approximation of smooth, 
non-analytic functions with Gaussians and other C°° positive definite functions poses a considerable 
challenge. In these cases, rates of decay are often known to be spectral [ID], but typically hold 
only for target functions that are infinitely smooth. More mainstream, linear error estimates for 
functions of finite smoothness (from Sobolev spaces, for example), where the L p norm of the error 
decays at a rate governed by the L p smoothness of the target functional have been more elusive. In 
this regard, the approximation power of the underlying spaces has only been thoroughly understood 
in the shift invariant setting, but that of the interpolants has not as yet been studied. In the case 
of Gaussians, the L2 error estimates fall under the shift invariant theory developed in (3j H] and 
generalized to L p by Johnson [7] Section 4]. 

1.1 Overview 

In this article, we demonstrate convergence rates for cardinal Gaussian interpolation for target 
functions having finite smoothness. The basic strategy we adopt has been developed in [13] and is 
based around a K- functional argument, usually brought about by band limiting the target function 
in a precise way. The techniques we use for estimating the error involve showing that interpolating 
by band- limited functions delivers precise approximation rates, and that such band- limited inter- 
polants form a very useful class of target functions on which cardinal Gaussian interpolation is very 
stable. In this case, the main challenge is to demonstrate this extra stability; this is accomplished 
by carefully controlling the multiplier norm of the Lagrange function. 

In Section [2j we describe in detail our main result and the strategy that we will use to obtain it. 
At the end of that section, we discuss a way to generalize the results to other, analogous situations. 
The approximation results for interpolation by band-limited functions is the focus of Section [3j 
The main tools for cardinal Gaussian interpolation and the key multiplier estimate are given in 
Section HI The extra stability results are demonstrated in Section [5j 

1.2 Notation and Background 

The symbol C, often with a subscript, will always represent a constant. The subscript is used to 
indicate dependence on various parameters. The value of C may change, sometimes within the 
same line. 

Let <S denote the space of Schwartz functions on IR n . The n-dimensional Fourier transform is given 
by f(0 = / K n f(x)e~ l ^' x > dx, and its inverse is / v (x) = (27r)~ n J Rn f(£)e l( ~ x '& d£. An important 
property of the Gaussian functions 

g : x h-> exp[— |x| 2 ] , (1.1) 

1 decaying exponentially fast 

2 This stands in contrast to nonlinear approximation, where the L p norm of the error decays at a rate dependent 
on a different smoothness norm, but generally Gaussians at different scales must be employed, cf. [6l [8] 
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is that they satisfy g = W 2 g(-/2). 

Given a multiplier m : M 71 — > C, which maps / to ^/ • m") the L p — > L p operator norm is denoted, 

ll m ll^1p := sup 
II/IIp=i 

Let Q C M n be a domain that satisfies a uniform cone condition. The Sobolev space Wp (ft) is 
endowed its usual seminorm | • an d norm || • defined by 

\f\w£(n) ■= sup \\D a f\\ Lp(n) , and \\f\\ W k {Q) := \\f\\ Lp (n) + \f\wp(ay 

\a\=k 

When just the seminorm is used, the resulting space is a Beppo-Levi space. We will denote it by 
Wp(fl). Finally, most of the time we will be dealing with ft = M. n . When that is the case, we will 
just use L p , Wp, or Wj; . 

We denote by B(x,r) the ball in W 1 having center x and radius r. The space of entire functions 
of exponential type, viewed as tempered distributions whose Fourier transform is supported in 
B(0, R), is given by 

PW(R) := {feS'\ supp(/) C B(0,R)}. 

This is the Paley-Wiener space of band-limited entire functions. In addition, we let PW£(i?) : = 
Wp n PW(-R), the band-limited entire functions in . 

2 The Main Result 

We consider interpolation of a continuous function / : C(W l ) — > C at gridded centers KIT 1 using 
elements of the linear span of shifts of a fixed Gaussian kernel g(x) = exp(— |x| 2 ), the span be- 
ing closed in the topology of uniform convergence on compact sets. In other words, we consider 
interpolation by functions of the form Sf t h(x) = ^2^ eh in a^g(x — (). 

The existence and uniqueness of the interpolant is a consequence of the existence/uniqueness of 
the Lagrange function Xh- Indeed, 

I h f{x):= J2 f(OXh(x-0- 

The Lagrange function is the function in the (extended) span of shifts of g equaling 1 at the origin 
and at all other dilated multi-integers. These have been studied in |14} I15| . as the cardinal 
interpolant: l}™\y) = J2jez n c j exp( — A|j/ — j| 2 ). The relation between the two is 

Xh ( X )= y, - o • ^ 

The problem that we have set forth for ourselves is to obtain a good understanding of how well 
Gaussian interpolants approximate functions in various Sobolev spaces. Our main result is the 
following: 
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Theorem 2.1 (Main Result). Let 1 < p < oo and k > n/p. There exists a constant C p so that 
for f € Wp the Gaussian interpolant I^f = X^e/iZ" f(OXh{' ~ satisfies 

\\f-hf\\ P <c p h k \\f\\ Wpk . 

For p = 1 and k > n or p = oo and A; > ; there is a constant C so that for f € W^iW 1 ) 

\\f -hf\\ p <c(i + \\ogh\Yh k \\f\\ w u. 

The strategy for proving the main theorem involves two steps, which we will discuss before we give 
the proof. The first is showing that interpolation of functions in W k by band-limited functions 
on KL n is possible, and that the band-limited interpolants approximate functions in W k very well. 
This is the content of the Approximation Property described in Lemma 12.21 below. We will prove 
this in Section El 

Lemma 2.2 (Approximation Property). Let b = it + e with < e < it. If f S W k for k > n/p, 
1 < p < oo then given h > there is a function f satisfying 

/€PW(6/fc) (2.2) 

such that 

f\hZ n = f\hZ n (2-3) 

\\f-f\\ Lp <Ch k \f\ Wpk (2.4) 

\f\ W H < C\f\ Wpk (2.5) 

The second step is to show that Gaussian interpolation is stable on the space PWj(6/fc), endowed 
with the Sobolev norm. We will prove this in Section [5j 

Stable Interpolation Property. Let k E N, k > n/p. We say that the interpolation operators 
Ih satisfy the stable interpolation property on the family PW k (b/h) if there is Q p : h — > (0, oo) so 
that one has 

\hf\w* < Qp(h)\\f\\w* for all f G PW k p (b/h). 
Proposition 2.3. Let 1 < p < oo. If the interpolation operators Ih satisfy the Stable Interpolation 

114/ -f\\ P < C(l + Q P (h))h k \\f\\ wi , / g W k . 



Property on PW k (b/h), then for k > n/p 



Proof. We note that, by the Stable Interpolation Property, ||/ — Ihf\\ P < ||/ — f\\ P + ||/ — Ihf\ 
since Ihf = Ihf '■ By the Approximation Property, the error ||/ — f\\ p is controlled by h k . Thus, 
the analysis of interpolation error / — I^f reduces to investigating its behavior on /iZ ra , where the 
error vanishes. An estimate on the size of a smooth function having many zeroes was proved by 
Madych & Potter Corollary 1]. Employing it, we obtain 



11/ - hf\\ P < Ch K \f - I h f\ ws < Ch« |/U + \Iff\wi 
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Thus, the interpolation error is controlled entirely by the norm of Ih as an operator from Wp 
Wp for k > n/p. Invoking the stable interpolation property, we obtain 



to 



\\f-hf\\ P < 



\\f-f\\ P + \\f-hf\\ P 

h k \\f\\ W k + h k (\\f\\ w u + Q p (h)\\f\\ Wpk 

h k {l + Q P (h))\\f\\ Wpk . 



< 



) 



< 



□ 



Proof of Theorem \2.1\ Finishing the proof only requires only showing that the Stable Interpolation 
Property holds, with the appropriate Q p (h). We do this in Lemma 15.41 There we also show that 



Apart from C p being explicitly dependent on p, the two constants depend on n, k, and the choice 
of the parameter b. Using this estimate on Q p in Proposition 12.31 then yields the result. 



Generalizations The interpolation problems considered here specifically involve only spaces of 
band-limited functions and spaces of Gaussians. However, it is worthwhile to broaden the context 
and describe these problems in a more general way. 

Let H C SI be a discrete set of points, which we will call nodes or centers, which will play the role 
of KL n above as sites for interpolation. Since S doesn't have to lie on a grid, we will describe how 
dense H is in SI using the fill distance or mesh norm, which is defined by h(E) := dist(x, S). 

Normally, one would not be dealing with all possible S, but rather with a specific class of sets, X. 

Suppose that for each H in X there are two spaces of functions, F3 and Ge- These form families 
.Fx = {Fk I £ € X} and Gx = {Ge I £ £ X}. Here, F^ and Ge are analogous to the band- 
limited functions PW (b/h) and the Gaussians, respectively. Both are contained in C(Q) n W^(O). 
In addition, we will assume that for each Ge there is an interpolation operator Ig B : C(O) — > Ge- 
Finally, we suppose that .Fx obeys an Approximation Property and Gx obeys a Stable Interpolation 
Property with respect to the family Fx- Then a nearly identical proof to the one for Proposition 
12.31 will establish this generalized version of that proposition. 

Proposition 2.4. Let Q C M. n be a region satisfying a uniform cone condition. For 1 < p < oo, 
let X be a collection of discrete subsets S C SI and let Fx be a family satisfying the Approximation 
Property. If the interpolation operators Ig B : C(ft) — > F"h satisfy the Stable Interpolation Property 
on Fx> then for k > n/p, 




1 < p < CO 
p = 1, oo. 



□ 



\\Ig-J - < C(l + Q p (E))h k \\f\\ W k , f e W^Q). 
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3 Band-limited Interpolation and the Approximation Property 



In this section, we will prove Lemma [2. 2\ which asserts that interpolation by band- limited functions 
satisfies the Approximation Property. 

Proof. Our aim is prove the existence of the interpolant /, along with the other properties. Let 
€ C°° be such that 

• supp C {£: |£|oc<7r + e} 
. 0(0 = 1 if |e|oo < 7T - e 

• £ £(£-27T/)=l. 



Note that if <J?(£) satisfies the first two conditions, then setting = fflf^jrj) defines a function 

jez 

which satisfies all three. Our candidate for / is the function 



.r 



9{x) := E ffafif 

jez n 

Clearly g is in PW(b/h) and also g\hz n = /Uz™, thus (|2.2p and (|2.3j) are valid for f = g. To see 
(j23j) and ([23]) also hold, write 

/ = /o + /i where /„(£) := /(0£(2/tf)- 
Define g = go + gi accordingly, i.e., 



and gi := g — go- Note that (for sufficiently small e) 



so 

Sd(ar) = /o0*0- (3-1) 
The Sobolev seminorms of go and /o are controlled by the seminorm of /: 

|/olw* < SU P P a (/*P)IU P < IIpIIlJ/Iw*, 1 <P<oo (3.2) 

|a|=fc 

where p(x) = m^yi 4>(x/2h) and \\p\\ii = \\</>\\ia is independent of h. On the other hand, since 
fx = f — / * p, we have, 

||/l|Up = ||/-(/*p)ilp<C^|/U, (3.3) 
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the usual error from mollification by a band-limited mollifier. What remains is to control the 
seminorm of g\ and to bound the error — /i|L- 

For |a| = r with < r < k, we have 

D a gi(x) = h(hj)D c 



When p = oo, it follows that 



|£> a <7l||oo < h~ r sup \h(hj)\ x sup V \[D C 



<Ch- r sup \fi(hj)\, 



j(rjri r- ' " ' * 1 

since the fact that is a Schwartz function implies that sup^gn X^ez™ l[^ a< ^] (v ~ j)\ — C- 
For p < oo, we have 



\D a gi(x)f = (h- r y\Y t fl(.hj)[D c 



< h- rp i 



1 + 



2n 



[D c 



-2np' 



P/P' 



Since is a Schwartz function, II (l + I ? — j\) 2n [D c 

comes about from the substitution y = ^ — j in the integration over M n . Consequently. 



f — •?') Hlp(r«) * s hounded by C r h n where h r ' 



\D a 9l \\ LP <h 



-^{^lA^l^l^xmax 



X 



-2np' 



P/P' 



-2n 



The sequence j i— >• (1 + ( y — j ' 1 is bounded in £ p , uniformly for all y G R n , by a constant C. 
Thus, for 1 < p < oo, 



Iflllwj < C r h r /i n ^ p ||/iUz' i ||£p(Z' 1 ) 
<c r ^ r fllMlLp + ^l/llw* 



(3.4) 



This follows by scaling the estimate \\F\zn \\e p (Z") < C(ll-^1lL p (R n ) + l-^1vK fc (R™))> which holds for 
k > n/p (this is a simple consequence of the Sobolev embedding theorem). By applying (|3.2j) and 
(13.31). we obtain 



\9i\w;<Ch k - r \f\ w ,. 



Finally, we have 



\9\w>f - \9o + 9i\wS - \ fo + 9i\wt - \ fo\wi; + bilwfc < C\f\ W k 
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which follows from (|3.2j) and above. So (|2.5p is established for 1 < p < oo. Also from (|3.ip . (|3.2p 
and (|3.4p we get 

11/ - 9\\lp = ||/i - < ||/i||lp + Hffill-W 
< Ch k \D k f\ LP . 

and (j2.4j) holds for 1 < p < oo. Hence, it follows that g is indeed /. □ 



4 Cardinal Interpolation with Gaussians 

We now investigate the Fourier transform of the Lagrange function, and, through this, the Fourier 
transform of the gridded Gaussian interpolant. Since XhiO = h n l)^(h ■ £) with A = h 2 , and 

exp(-|e| 2 /(4A)) 



because 



we have 



XhiO 



£ fee z«exp(-(|£-27rA;|)7(4A)) 
expH|£| 2 ) 



£ fce ^exp(-i|£-^| 2 ) 



(4.1) 



Throughout the rest of this article, the multiplier is the subject of much of our investigation. 
Its multiplier norm controls the the stability of the interpolation process on Sobolev spaces, and 
ultimately, it determines the rate of decay of interpolation error. Because it is, evidently, a re- 
fold tensor product of univariate multipliers: rn)?\£) = mj^(£i) • • • mj^(£ n ), many of our results 

can be obtained by considering the 1 dimensional case, where we suppress the dimension to write 

[i] 

m h := m l h '. 

Finally, formula (j4. 1 j) allows us to easily express the Fourier transform of the interpolant: 



WiO = \52f(hj)xh(--hjj\ iO 



4%) 



/3G27rZ 



mt\0- 



(4.2) 



(The final equality follows from Poisson's summation formula.) 



4.1 The univariate multiplier ran 



In the interest of keeping results self contained, we now provide some simple estimates on the size 
of the multiplier. Such estimates (and stronger ones) could, with some effort, be gleaned from the 
work of Riemenschneider and Sivakumar. However, the ones we provide here are easier to obtain 
than those presented in [151 Theorem 3.3 and Corollary 3.4], yet totally sufficient for our purposes. 
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We proceed in two stages. In the first stage we obtain estimates on that hold independently of 
h, but are rather slow (in particular, they do not show that m/ l v is integrable). In the second stage, 
we demonstrate that m/j V (£) decays like 0(|x| -2 ), and is, hence, integrable, but these estimates 
depend strongly on h. 



First estimate of m/j V : Since m/ l (^) > 0, we have |m/ l v (x)| < m/ l v (0) = 1. On the other hand, 
from (14. ip . 



-i 



-i 



l + 2^e"^"cosh(vr^//i) =: (d (0) 



k=0 



Hence m' h (£) = -(m fe (^)) 2 di(^) and m{£(£) = 2(m /l (^)) 3 (di(^)) 2 - (m fe (£)) 2 n 2 (£), where we have 
defined c£i(£) := d' (£) and ^ := n[(^). Hence, 



fc=l 



irk 



irk£ 



oo / , \ 2 



E/ 71 ft; \ 7t k 
I — I e ft 2 cosh( 



fe=i 



It follows that < for positive £, and, by symmetry, m' h (£) > for negative £. Therefore, 

Jo°° Kh(f )l = - /o°° d £ = "»ft(0) < 1. It follows that 



\m' h \\i < 2, 



which implies that |m/ l v (x)| < 2/|x[. Thus, 



\ m h V ( x )\ < min(l,2/|x| 



(4.3) 
(4.4) 



Second estimate of m/ l v : Using some simple algebraic manipulations, we may rewrite the second 
derivative as 



m h{i) = m h(0 



Z^fcez I h 



2 lie 27rfc|2' 

I e 4l« hi 



v-^ _I|tf_2£fc|2 

E fce z e 4 ^ h 



e 4 



(4.5) 



This leads us to the the following estimate: 



Lemma 4.1. Let /i < 1. There exists a constant C so that for k = 1,2, ... , and £ € [ t , 



K(0I<C U ro fc (fl. 



Proof. We write 



z^fcgzl— J e 4 " fe 



7 := 



. Efc G Z e 4 ^ h 



and II : = 



We split the numerator of I to isolate its two principal terms 



fcez 



/C7T 



lie 27rfc |2 



< 



kn 



.i |c _M|zll |2 



+ e 4 



I |f 27rfc |2 



fe=l 



/C7T klT \ \irk\2 

~~h + ~h) ' 
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Since it is a series of nonnegative terms, the denominator can be bounded below by its two largest 

1 If 27rfc |2 l i- 2_7r(fc— 1) 2 1 i f 2-n-fc |2 71-2 

summands: 2^ g ^e 4 |? h l > e * IS & +e h I > 2e ih 2 ". Therefore, 




/ < 

A similar argument shows that 

( kn\ ( kn kn\ _( , |»*|2_ _zdM 



fe=i 



□ 



An immediate consequence is that for h < 1, < Ch 3 . Indeed, one may estimate the 

integral on an interval around the origin [— 2ir/h, 2tt, h] and along the punctured real line f2 := 
E \ [-2vr//i, 2vr//i] to obtain 

r-2n/h 

|m"(0|d£ < C/h 3 

-2n/h 

and 

It follows that 

/i 3 |x| 2 



m, v (x)| < T ^ i? . (4.6) 



4.2 The multiplier norm of mj^ 

In the Section [5] we investigate the Sobolev stability of Gaussian interpolation on spaces of band- 
limited functions. Of particular importance are the operator norms of the convolution operators 

[nl 

with Fourier multiplier m h . These can be estimated in the case p = 1 and oo by using the bounds 
on the cardinal interpolant obtained in |15} Section 3] (although the ones developed in the previous 
subsection, (|4.4p and ()4.6f) . suffice). In case 1 < p < oo, multiplier norms are estimated by appealing 
directly to the formula (14. ip . and making a comparison to the Hardy-Littlewood maximal function 
and the (maximal) Hilbert transform (this is a continuous version of an idea used in [121 Theorem 
3.1]). 

Lemma 4.2. For 1 < p < oo, there is a constant C p so that 

\\m [ Z ] \\ Mp <C p . 
For p = l,oo there is a constant C so that 

irfll^i = <C(1 + | log /*[)». 
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Proof. Because mi" is a tensor product of univariate multipliers, it suffices to consider the case 
n = 1. 

For l/p+ 1/p' = 1, we have \\m h \\ Mp = \\rn h \\ Mpl . Thus, \\m h \\ Ml = \\m h \\ Moo < ||m h v j|i, by 
Holder's inequality. ^From (|4.4p and (|4.6p . we have 



|m/ l v (x)| dx < 2 



1 + 2 / (x) _1 dx + Oh' 1 / (/ix)~ 2 dx 



[C + 12\\ogh\] < C(l + \logh\) 



We now consider the case 1 < p < oo. Let / £ L p and let g € L p / with ||/|| p = 1 = ||g|| p '. We 
can estimate ||m/i||x p by the supremum of the expression | f R f R f(x)mf^ ' (t — x)g(t) dx dt . To this 
end, let n h (t) := R\ (t- h,t + h). Then, 



f(x)m h y (t - x)g(t) dx dt 



< 



t+h 



f(x)m h y (t - x)g(t) dx dt 



t-h 



+ 

' iRJn h (t) 
--: I + 11 



f(x)m h v (t - x)g(t) dx dt 



The first integral I < J R h~ l f*Z£ \f(x)\\g(t)\dxdt, can be compared to an integral involving the 
maximal function of /, /"(i) := sup e>0 (2e) _1 \f(x)\ dx. Thus, 

I<2 [ f$(t)\g(t)\dt<2\\f%\\g\\ p , < 2C p \\f\\ p \\g\\ p ,. 

JR 



To treat II, we note that, m' h (^) is integrable, and II can be estimated by employing the Fourier 
transform of rrih V , (this is a trick similar to the one used in |12j): 



II 



uJa h (t) 



f(x)m h (Oe i{t - x)i g(t)d^dxdt 



_I^Lm' h (Oe< t - x) h(t)d^dxdt 
m Jn h (t) Jr H x _ *J 



/(*) 



n h (t) - 1) 



^ l - x ^g{t)dxdtd^. 



The second equality follows by integration by parts, while the third follows by Fubini's theorem, 
since m' h is integrable on K and S^L g(t)lci t (x) is integrable on M 2 . It follows that 



II < 



"4(01 



n h (t) 



e-^ x f(x) 
x — t 



dt 



\g(t)\dtd£ 



< [ \m' h (0\ [ [H(e-*/)](t) \g(t)\dtdC 

JR JR 

< c p\\m' h \\i\\f\\p\\g\\ p > <2C p \\f\\p\\g\\p>. 
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In the second inequality, we use the maximal Hilbert transform [HF] (x) := sup e>0 | J n n\ F(t)-0^\. 
The third inequality follows from the fact that the maximal Hilbert transform is strong type (p, p) 
(i.e., ||HLF||p < Cp||F|| p ) for 1 < p < oo, as observed in [1, Theorem 4.9]. The fourth inequality 
follows from (JOJ). □ 



5 Stable interpolation of band-limited functions by Gaussians 

We consider in this section interpolation of functions in PW(6//i) and show that Gaussian interpo- 
lation restricted to such functions is stable with respect to each Sobolev norm Wjj, with k > n/p. 

Before stating and proving the result in its full generality, we indicate how the proof works in the 
univariate case. We use (|4.2p to write, for a nonnegative integer k, 



(4/)°'(f) = e 



/3e2nZ 



(0 = £ G iA0 



k 



where G jiP := LP [(e^W^f) * (ro A ) v ] . Clearly, < Ej=oE,3 \\GjAv 

The result we are after requires estimating ||G-j a[L for various values of j3. These estimates fall 
into three categories: /3 = 0, |/3| = 27T, and |/3| > 2ir. 

Estimating ||Gvj,o||p : In this case, * (m/ l ) v || < ||-D J '/llpll m /i||.Mp from which we obtain 

T^j<k \\ G jfl\\p - C \\f\\w k \\ m h\\M P - Thus, from Lemma[01 



Eii^ii p <{^ 



'pWIWw* 1 < p < 00 , g 1 s 

+ |log^l)ll/llw* P=l,oo 



J 

where C and C p depend only on k and p. 

For some of the other terms, we need to use a special cutoff function. We first consider a smooth 
bump function, v. It is a non-negative C°° function satisfying 

• v(£) = if |f | > 2e and u(f) = 1 if |f | < e, 

• v(-£) = v(£). 

We use v to construct univariate cutoff functions cp having support in [— tt — 2e,ir + 2e]: 

• (p(t) = 1 for — 7r < t < it; 

• <p(—t) = ip(t) = v(t — tt) for t > it. 

Estimating HGj^Hp, for j/3| > 2ir: In this case, the fact that / is band limited, supp/ C 
-6(0, (tt + e)/h), allows us to write 
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The norm of the multiplier t\ := Tij j( g(£) := £V(M£ — T;)) m h(Q can be estimated by ||ti||i, which 
we can estimate using the fact that 



max 



^(0 



dC / \g(0\ <%) <k 

Jm J 



\g{x)\ < 



K 



(1 + M) 2 

By (|5.9|) of Lemma 15,51 we have the following. 

Claim 5.1. For f3 G 2xZ n , \/3\ > 2tt, \\T X j,p\\i < (C//i)|/3//i| (2+j) exp(-c|/3| 2 //i 2 ). 



\g\\l < CK. 



Therefore, 



IGj^Hp < ll r i||i 



h 



< 



(2+i) 



exp 



h 2 



(5.2) 



Estimating \\Gj b\\ p , for |/3| = 2tt: This remaining case is quite similar to the previous one. We 
again use the fact that / is band limited, although we need to exercise extra caution since the 
cutoff <p(h(£ — (3/h)) overlaps a narrow region (near to /3/2h) where mh(£) is close to 1 and is 
very large. 

When (5 = 2ir we write ip(t) = u:(t) + v(t + ir), with oj having support on the (non-symmetric 
interval) [— 7r + e, 7r + 2e] and equaling 1 on [— 7r + 2e, 7r + e]. (When (3 = —2tt an obvious modification 
cp(t) = Cb{t) + v(t + 7r) is made.) This allows us to write 

= f^-^)ev{h(C-^))m h (0 + f(C-l)e^-27T/h))m h (0 



-■ ft-^(0 + /({-^(f). 



(5.3) 



We first investigate T2 = T2j t p by rewriting the monomial £ J as a Taylor series about 2ir/h, obtaining 
r 2 (0 = ev(h(C - l))m h (0 = £ ( J p ) (t - V(£ - l))m h ® 

x m h {Q x 



The multiplier norm of /i® := ELo © (x)'"' = ELo ffl (which is 

a Schwarz function, since the support of v(-— 7r) is positive distance from 27r) is uniformly bounded 
(in /i) by the constant 

: = E G) ( 2 ^'~ 

£=0 v 7 

Thus we have shown the following: 

Claim 5.2. The multiplier T2j p(£) can be written T2j = (£ — fi/h) 3 x mh{i) x where /x 

has multiplier norm \\fJ-\\M P < bounded independent of h. 



v • — 7T 
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The multiplier t^j p is controlled in a similar way to t\. The only modification is that we use 
estimate (|5.1ip of Lemma 15.51 utilizing the fact that uj{- — 2tt) has support a positive distance 
(namely e) from tt. In particular, it satisfies condition (|5.1U|) Thus, we obtain 



< 



i+3 



,-\c\ 2 /h 2 



and we have demonstrated the following claim. 

= 2ir. There is a constant C depending only on j for which the multiplier T3 



Claim 5.3. Let 

satisfies 



ll(T3 Jl/3 ) v ||i<c. 



For \/3\ = 2ir, applying the Claims [5721 and [5T51 to (|5.3p . it follows that 



\G jiP \\ p < \\(r 2 f(- - P/h)) v \\ p + \\T 3 f(- - P/h)) v \\ p 

v 



< lCj\\m h \\Mp 

< fcj\\mh\\M P 



p/hyf(--p/h) 
+ c\\f\\ p 



+ c 



f{--P/h) 



(5.4) 



Summing ||Grj j( g|| p over < j < k and (3 E 27rZ, and employing the estimates (15. ID . (15. 2|) and ([5] 
we observe that \\hf\\w£ < c (\ + ll m fc|I.M p )||/||w* for / G PW((7r + e)/h). 

We now give the general, multivariate result. 

Lemma 5.4. Let < e < ir/2. Cardinal interpolation by Gaussians on KL n restricted to band- 
limited functions in PW((7r + e)/h), satisfies the Stable Interpolation Property, with 

• Qp(h) < C p , a constant depending only on e, p, n and k when 1 < p < 00, 

• Qp{h) < C(l + I log/i|) n , a constant depending only on e, n and k when p = 1, oo ; 



In o£/ier words, 



\hf\wi < QpWII/llw* /° r / G P( 



7r + e. 



Proof. We use (|4.2p to write, for a multi-integer |a| < k, 



where G, 



a,p 



D c 



D a hf(0 



(feWih-)) * ( m N)Vj . Clearly, < £ H < fc Ep \\Ga,p\\p> and the remain- 

der of the section is concerned with estimating HGq,^ || p for various f3. 

We write Gp i0l (£) = M a (£) f (£ — f3 / h) employing the tensor product multiplier M a (£) := £ a rn^(£,). 
We estimate these terms by observing that the support of /(£ — /3/h), which is a neighborhood of 
/3/h, and therefore lies in a region where \i a is small (for most values of j3 7^ 0). 
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This heuristic is complicated by the fact that, for certain values of (3, the neighborhood of (3/h 
overlaps a region where mj^ is near to 1 and may be quite large. Therefore, we must be 
somewhat careful. 

Fix (3 = (Pi,... ,P n ) € 2nZ n . Before proceeding, we partition (1, . . . ,n) into three subsequences 
(Ji, . . . , J ni ), (Ki, . . . , K n2 ) and (L 1? . . . , L„ 3 ) , with n x + n 2 + n 3 = n, and where 

• (Ji, . . . , J ni ) are the indices J where \(3j\ > 2tt 

• (Ki,...,K n2 ) = supp(/3) , 

• (Li, . . . , L„ 3 ) are the indices where = 2ir. 

As an example, in dimension n = 5, we have for j3 = (— Air, — 2ir, 0, 0, 6ir) that (Ji, J2) = (1,5), 
(Ki,K 2 ) = (3,4) and L x = 2. 

Then, because the multiplier and the monomial £ a are tensor products and can be factored 
over the three subsequences we have just constructed, and because because / is band limited, we 
have 



M a (om-^) 



3=1 



{[(Uf'^'n^K 



Y\(^ L] r^(h^ Lj -^f))m h (i L 



x fit - f ) 



In other words, M a (£)/(£ — ^) can be written as a product of tensor product multipliers applied 
to /(£ - f ), namely 



M Q (0/(£-f) = M0]x 



TL2 



We have written ^tf) := < 7 1)a , j8 (0 := X^U^T^^K^ ~ jt))™^) = U/ .r-JU I In 
a similar way, we identify the factor where ft = as n™=i(£K j ) aKj mh^Kj) and the factor where 

n ^ 



\f}\ = 2n as a 2 (0 := <x w (0 := 1 1? ; : ^ /- r (/'(£/.. - ^m^)- It follows that 



IG 



a,pilp 



(M a f(--P/h)) 



(<T2,a,f>x{-) aK f(--P/hj) y 



(5.5) 



where we write (-) ax : £ ^ £ QK == Il?=i 
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Estimating \\(o- 1 ^ l3 ) y \\ Mp : We use the rough estimate \\(cri :Qil3 ) v \\ Mp < ||(o-i iQi/ g) v ||i. The Li 

norm of each univariate factor ti(£j.) := (£,j j ) ajj i p{h(^j j ^-))m/ l (^j j ) is bounded by Claim loTTj 

from which we obtain 



Estimating 



(<n) V Hi< 



<r W x(r/(-W 
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c_ 



n 

3=1 



fa 



h 



a j . +2 



(5.6) 



There is an immediate decomposition of <X2(£) into 



nweii)+7s(M- 



Thus Claims 15.21 and 15.31 and the fact that 
\ W k imply that 



< 



ff 2 x (.)"*/(■ -/3//i) 
with constant C depending on n and a. 



< C 



"3 



nt 1 + iwuo 



(5.7) 



Applying estimates (|5.6p and (|5.7p . which control ||(o- 1)Q , j( g) v ||_A^ p , and 
respectively, to (|5.5|) . we bound the sum of the ||G a)J g||p's: 



(<r W xQ a */(--/3/^)) 



\\G a> p\\ p <C(l + \\ mh\\M p ) 7 

/3e2vrZ™ 

and the Lemma follows from Lemma 14.21 



n s+ 2 E 

j=i \ 



e=2 



2tt£ 



OLA+2 „ X 
J 2 ^ I 2 

e c l h I 



(5.8) 



□ 

Lemma 5.5. Suppose that e > and i/iai (p is a C°° function with support in B(0,tt + e). If 
f3 € 27rZ, f3 > 2ir and < h < 1 f/ien t/iere exist constants c, C > 0, depending only on e and k so 
that 





d 2 


1 

m 






= 2tt 



£V(M£-f)K(£) 



^<ch- 1 ^ 2+fc exp 



supp(0(- - p)) n 5(0, vr + e) = 0, 
t/iere exist constants c, C > 0, depending only on e and k so that for < h < 1 







/ 


d£ 2 



d£ < C/i" (3+A;) exp 



(5.9) 



(5.10) 



(5.11) 
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Proof. Let b = it + 2e. We prove this by making the estimate 





d 2 - 


/ 


de 



«<r? 



max 



de 



By applying the product rule to the expression being maximized, we obtain 



de 



Y, c 7 x x - -))) x (^m^e)) 

l7l=2 ' V 7 



with C. 



2! 



7 71172173! ' 



Without loss, it suffices to consider only components where both derivatives are on the third factor, 
i.e., those of the form £, k (f)(h(t; — ^)) > since the other terms are small compared to these. 

Indeed, |£>Tf*| = C|£ fc -7| < C\f3/h\ k for £ G [(|/3| - 6)//i, (|/3 + 6)/fc] (with C depending only on fc), 

<HM£ - ?)) < C7*' 7 ', since T : = £)7 



and max, g+s i 
is a C°° function with the same support as 

Hence, 



max, r£-6 g+6 i 

? fc l h ' h J 



max 



de 2 



< c 



max 
heesupp(^(— y8)) 



The result now follows directly from Lemma 14.11 Indeed, we have: 



max 
fc£esupp(#-/3)) 



K(0| <c 



+ 6 



max mrihiCj < C 

h€68URpW(-^)) 



|/9[ + & 



e 4 



^l 2 



max 



hfGsupp(^(-/3)) 2e -i|f I 



If (|5.10p holds, then the expression being maximized can be controlled by exp(— j ^ n+e ^ ~ ^ and 
(gUD follows. 



On the other hand, if \/3\ > 4-7T, then /i£ 6 supp(^(- — /?)) implies that |£| 2 < 



exp(-i|£| 2 )exp 



1 7T 

4 Ul 



< exp(— c^Jr), from which (|5.9p follows. 



and 
□ 
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